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Abstract 

We give a new proof of correlation estimates for arbitrary moments of the resol- 
vent of random Schrodinger operators on the lattice that generalizes and extends 
the correlation estimate of Minami for the second moment. We apply this mo- 
ment bound to obtain a new re-level Wegner-type estimate that measures eigenvalue 
correlations through an upper bound on the probability that a local Hamiltonian 
has at least n eigenvalues in a given energy interval. Another consequence of the 
correlation estimates is that the results on the Poisson statistics of energy level 
spacing and the simplicity of the eigenvalues in the strong localization regime hold 
for a wide class of translation-invariant, selfadjoint, lattice operators with decaying 
off-diagonal terms and random potentials. 



1 Introduction: Correlation Estimates and Energy 
Level Statistics 

Correlations between various families of random variables associated with disordered sys- 
tems are an important aspect governing the transport properties of the system. For ex- 
ample, the conductivity is expressible in terms of the second moment of the one-electron 
spectral density. Another example is the correlation between the energy levels of non- 
interacting electrons for finite-volume systems and their behavior in the thermodynamic 
limit. Some of the first studies of energy level correlations were made by Molchanov [Hj 
and by Minami [13] for systems in the strong localization regime. Molchanov [H] studied 
a family of random Schrodinger operators in one-dimension with a random potential given 
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by q(t,u) = F(x t (uj)), where x t (uj) is Brownian motion on a compact manifold K and 
F is a smooth, real-valued, Morse function on K. It is known that this model exhibits 
Anderson localization at all energies (cf. [15111]). Minami [13] studied the lattice Ander- 
son model (see ([I])) in any dimension with a bounded random Anderson-type potential 
for energies in the strong localization regime. These authors proved that, under certain 
hypotheses, the normalized distribution of electron energy levels in the thermodynamic 
limit is Poissonian. This is interpreted to mean that there is no level repulsion (nor at- 
traction) between energy levels in the thermodynamic limit provided the energy lies in the 
strong localization region. This is in contrast to the expected behavior when the energy 
lies in the region of transport and strong correlations between energy levels are expected. 
In this case, the expected eigenvalue spacing distribution is a Wigner-Dyson distribution 

(cf. [E]). 

The precise formulation of this result is as follows. The standard Anderson model 
studied by Minami is given by the following random Hamiltonian acting on £ 2 (Z d ), 

H^(x) = 1>(V) + V(x)ip(x), xeZ d , ^ef(Z d ), (1) 

y;\y-x\=l 

where the potential to = (V(x)) xeZd is a family of independent, identically distributed 
random variables with common distribution with a density p(V(0)) such that ||p||oo = 
supy( ) p(V(0)) < oo. Let H\ denote the restriction of to a box AcZ d with Dirichlet 
boundary conditions. The spectrum of Ha is finite discrete and the eigenvalues E^iui) 
are random variables. For any subset J C R, we let Ea(J) be the spectral projection for 
Ha and J. The integrated density of states (IDS) N(E) is defined by 

|A|^oo |A| 

when this limit exists. It is known (cf. [3J [TB]) that, for the lattice models considered 
here, this function exists and is Lipschitz continuous (at least if the density has compact 
support). Consequently, it is almost everywhere differentiable, and its derivative n(E) is 
the density of states (DOS) at energy E. In order to describe energy level correlations, 
we focus on the spectrum near a fixed energy E. We define a point process £(A; E)(dx) 
by 

Z(A;E)(dx) = J2 S((\A\(Ef(u;)-E)-x)dx (3) 

jeN 

The rescaling by the volume |A| reflects the fact that the average eigenvalue spacing is 
proportional to |A| _1 . Minami proved the following theorem. 

Theorem 1 Consider the standard Anderson model (TJP and suppose the DOSn(E) exists 
at energy E and is positive. Suppose also that the expectation of some fractional moment 
of the finite-volume Green's function decays exponentially fast as described in l[21\) . Then, 



2 



the point process ([3J) converges weakly as |A| — > oo to the Poisson point process with 
intensity measure n(E) dx. 

Minami's result requires two technical hypotheses: 1) the density of states n(E) must 
be non-vanishing at the energy E considered, and 2) the expectation of some fractional 
moment of the finite- volume Green's function decays exponentially. Wegner[2"T] presented 
an argument for the nonvanishing of the DOS n(E) and a strictly positive lower bound 
was proved by Hislop and Miiller [8] under the assumption that the probability density 
satisfies p > p m i n > 0. Suppose the deterministic spectrum of H w is [£_,£+]. Then, for 
all e > 0, there is a constant C e > 0, depending on p m i n , such that n(E) > C e for all 
E G [XL + e, S + — e]. Exponential decay of fractional moments of Green's functions for 
random Schrodinger operators was established in certain energy regimes by Aizenman and 
Molchanov [2], by Aizenman pQ, and by Aizenman, Schenker, Friedrich and Hundertmark 

Additionally, Minami's proof rests on a certain correlation estimate for the second 
moment of the resolvent. It is this estimate that interests us here. We present a new 
proof of this estimate that generalizes Minami's result in two ways: 1) it holds for gen- 
eral bounded, selfadjoint Hamiltonians Hq, including magnetic Schrodinger operators and 
operators with decaying, off-diagonal matrix elements, and 2) it holds for higher-order 
moments of the Green's function. These generalizations of Minami's estimate were re- 
cently also obtained by Graf and Vaghi [7] with a different method, which we outline in 
section 14.21 below. 

We also apply this moment bound to prove a new estimate on the probability that there 
are at least n eigenvalues of a local Hamiltonian in a given energy interval. We interpret 
this as an n-level Wegner-type estimate that bounds the probability of n eigenvalues being 
in the same energy interval. As such, it is a measure of the correlation between multiple 
eigenvalues. 

Minami's estimate may be stated in several ways. For z G C + , we let R A {z) = (H A — 
z)^ 1 denote the resolvent of the finite- volume Hamiltonian on £ 2 (A). The corresponding 
Green's function is denoted by G A (x, y; z), for x, y G A. Minami stated the estimate this 
way (Lemma 2, [13J). 

Lemma 1 For any z G C + , any cube A C Z d , and for any x, y G A with x ^ y, we have 



E 



det 



3Ga(i,i; z) ^sGa(x, y; z) 
%G A (y,x;z) %G A (y,y;z) 



<^\\P\\L- (4) 



However, for many purposes, it is clearer to note that terms as on the left side of (j4]) 
arise when evaluating (Tr{QRA_}) 2 — Tr{(QR\) 2 } in the canonical basis of £ 2 (A). Thus 
(j3J produces the bound 



E 



(Tr{QR A (z)}y-Tr{(ZR A (z)y} < ^||p||^|A|^. (5) 
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(6) 



for any interval J C R. This estimate ([HD was used by Klein and Molchanov [TU] to 
provide a new proof of the simplicity of eigenvalues for random Schrodinger operators on 
the lattice, previously shown by Simon |20j with other methods. In fact, from Chebyshev's 
inequality, we can write (jHJ) as 



It is crucial for the applications that in the bound of the left side of (j7|) the exponents 
of both the volume factor and the length of the interval \ J\ be greater than one. We 
mention that a bound of the type ((Tj) is not known for random Schrodinger operators 
on L 2 (A), for A C M d . This is the main remaining obstacle to extending Minami's 
result on Poisson statistics for energy level spacings, and the Klein-Molchanov proof of 
the simplicity of eigenvalues, for energies in the strong localization regime, to continuum 
Anderson- type models. Our original motivation for this work was two-fold: First, we 
wanted to find another proof of Minami's miracle in order to better understand it, and, 
secondly, we wanted to try to generalize the Minami estimate so that it was applicable to 
other models. 

1.1 Contents of this Article 

We state our main results, a generalization of Minami's correlation estimate, Theorem [21 
and its application to an n-level Wegner estimate, Theorem[3l in section 2. Given Theorem 
[2j we prove the n-level Wegner estimate in section 2. We prove the main correlation 
estimate in section 3 using a Gaussian integral representation of the determinant. In 
section 4, we discuss applications to energy level statistics and the simplicity of eigenvalues 
in the localization regime for general Hamiltonians, and a related proof of the correlation 
estimate due to Graf and Vaghi [7]. We also discuss the related works by Nakano [15] and 
Killip and Nakano [9] on joint energy-space distributions. For convenience, a proof of the 
Schur complement formula is presented in the appendix. 



We consider random perturbations of a fixed, bounded, selfadjoint, background operator 
Hq. The general Anderson model is given by the following random Hamiltonian acting 



F{TrE A (J) > 2} < y ||p||LM 2 |A| 2 . 



(7) 
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2 The Main Results 



on £ 2 (Z d ), 
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H u if)(x) = H ip{x) + V{x)ip(x), xeZ d , ^ef(Z d ), (9) 

that generalizes (JT]). The potential uj = (V(x)) x&zd is a family of independent, identically 
distributed random variables V(x) with distribution given by a density p(V(0)) such that 
IIpIIoo = sup v(0) p(V(0)) < OO. 

Among the general operators Hq, we note the following important examples. The first 
family of examples include nonrandom perturbations of the lattice Laplacian L, defined 
by, 

{I4)(?)= E 1>(V), (10) 

y;\y-x\=l 

by T C Z d -periodic potentials V on Z d so that H = L + V . Here, the group T is some 
nondegenerate subgroup of Z d like NZ* d , for some N > 1. The second family of examples 
consists of bounded, selfadjoint, operators H with exponentially-decaying, off-diagonal 
matrix elements. The third family of examples are discrete Schrodinger operators with 
magnetic fields, 

{H^){x) = £ (^)-e^V(y)), (11) 

y;\y-x\=l 

where A(x,y) = —A(y,x) is nonvanishing for \x — y\ — 1 and takes values on the torus. 
The operator Ho need not be a Schrodinger operator but simply a bounded selfadjoint 
operator for Theorems (2J [3], and HI The boundedness of H is not essential, but we will 
require this in order to avoid selfadjointness problems. 

When we consider localization and eigenvalue level spacing statistics in section 4, we 
will require that, in addition, the selfadjoint operator Hq is translation invariant with 
the off-diagonal matrix elements, |(x|f/o|y)|, decaying sufficiently fast in \x — y\. We will 
discuss the required properties of Hq further in section HI 

2.1 Generalization of Minami's Correlation Estimate 

For any subset A C Z d , we define Pa to be the orthogonal projection onto £ 2 (A), so that 
P\f(x) = J2 y eA f(y)$x,y By H\ we denote the restriction PaH^Pa of to £ 2 (A). 
Let A C A and note that V commutes with Pa'- Va = PaV = VPa- For z G C, with 
$S(z) > 0, the matrix-valued function Qa{z) = Pa(H\ — z)~ 1 Pa has the following property 
(see [13J for the case n — 2). 

Theorem 2 For Q(z) > and any subset Ac A, with |A| = n, the following inequality 
holds: 

E(det{% A (^)})<7rlp||^, |A|=n. (12) 

A new proof of this result, using the representation of the square root of a determinant 
by a Gaussian integral, will be given in section [3j It is a generalization of Minami's result 
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Lemma [I] where H is the discrete Laplacian L, defined in ( TTUj) . and n = 2. In the case 
n = 2, we may write Pa = \ x )( x \ + \y)(y\i f° r ^ Vi so that 

„ ^ _ f G A (x,x;2;) G A (x,y;z) \ . , 

^ z >-{ G A (y,x;z) G A (y,y;z) ) ^ 

where, as above, G\(x,y;z) is the Green's function for H\. Thus Lemma [1] follows 
from Theorem [2] if we note that in this case one has G A (x,y;z) = G A (y,x;z) and thus 
^9 a = (flA — 9a)I^ 1 m CGD is the same as the matrix on the left of (jlj). 



2.2 The n-level Wegner Estimate 

We use Theorem [2] to prove a new estimate about multiple eigenvalue correlations. We 
begin with the observation that the left hand side of ffT2l) can be interpreted in terms 
of the eigenvalues of an operator acting on a certain antisymmetric subspace of a finite 
tensor product. 

Theorem 3 Let A = A* be a self adjoint operator on £ 2 (A), A C Z d finite, with eigenval- 
ues a\ < a2 < • ■ ■ < on, where N = |A|. Then, the following holds 

^2 det{P A AP A } = ^ a h ...a in . 

AcA;|A|=n l<h<-<in<N 

Moreover, if7i n = £ 2 (A) An is the n-fermion subspace, let A An be the restriction of A® n to 
H n . Then 

det{P A AP A } = Ti Hn (A An ) . 

AcA;|A|=n 

Proof: The first identity is a trivial consequence of the second. For indeed, the eigen- 
values of A An are products of the form . . . ai n with 1 < i\ < ■ ■ ■ < i n < N and the 
trace is the sum of the eigenvalues. 

To prove the second identity, for each x G A, let e x be the unit vector in H\ = £ 2 {A) 
supported by x, namely e x (y) = S XjV . Then {e x ; x G A} is an orthonormal basis of 
TCi- Let A be ordered so that we may write X\ < x 2 < ■ ■ ■ < X]y. Then, {e Xii A ■ ■ • A 
&x in I Xu G A, 1 < ij < N} is an orthonormal basis of Ti, n if we restrict to indices so that 
Xi x < Xi 2 < ■ ■ ■ < Xi n , with 1 < ij < N. Thus, the trace on TL n can be expanded as 

Tr Wn (A An ) = ^ A ■ ' ' A e ^\ AAn e Xl A-AeJ. 

Xl<--<X„ 

If A = {xx, ■ ■ ■ ,x n } (where the labeling is such that X\ < x 2 < ■ ■ ■ < x n ), then the 
definition of the determinant gives 

(e Xl A • • • A e x JA An e xl A • • • A e Xn ) = (e Xl A • • • A e Xn \(P A AP A ) An e Xl A • • • A e Xn ) 

= det{P A AP A }. 
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□ 

We can now combine the previous two Theorems to generalize (GO) and prove the 
following n-level Wegner estimate. We point out that this estimate holds for all energy 
intervals J in the spectrum of H A . 

Theorem 4 For any positive integer n, and interval J C K, and any cube A C Z d we 
have 

F(TrE A (J)>n)<^\\p\UJr\Ar. (14) 

Proof: By taking A = QR A (z), Qz > 0, in Theorem [3] and using the result from Theo- 
rem [2] we get 

E(Tr Hn ((ZR A (z)) An )) = Yl E(det{% A (^)}) 

AcA; |A|=n 

< ( |A 'Vl*- (15) 



n 



For ( = a + ir, a G R, r > 0, define 

/<(*) 



(x — a) 2 + r 2 



If (a,b) C J C [a, b] and z := (a + b + i\J\)/2, then Xj{ x ) < \Afz{x) for all x G R and 
thus, by the spectral theorem, 

E A (J) < \J\QR A (z). 

This carries over to Ti, n as 

E A (J) An < \J\ n (SR A (z)) An . (16) 

If X is the range of E A (J), then E A (J) An is the orthogonal projection onto the subspace 
X An of H n . Thus 

AlJJ j | Q if TrE A (J) < n. 

Chebyshev's inequality and the elementary fact that k/n < for all > n imply that 

P(Tr£ A (J)>n) < ((Tr£ A ( J)) ■ X {^E A (j)>n}) 

( {TrE A (J)\ 
~ V V n / X{Trf?A(J) - n} 
= E(Tr Hn E A (J) An ). 
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Using the bounds ATE]) and (fi"5l) finally yields the desired result, 
P(Tr£ A (J) > n) < \J\ n E(Tr Hn {QR A (z))' 



< \J\ n ( lAl U n \\p 



n 



"■II 

oo 



7T™ 

< \\n\\ n / " A " 



□ 



3 The Generalized Minami Correlation Estimate 

Our approach to the Minami correlation estimate Lemma [U and its generalization, is to 
work with the resolvent, rather than the Green's function. We use the Schur complement 
formula to isolate the random variables and the representation of the inverse of the square 
root of a determinant by a Gaussian integral, see (fT8j) . The proof of Theorem [2] requires 
several steps. As in section |2j we let A C A and denote the orthonormal projection onto 
£ 2 (A) by P A - We define H A = H A - V A , and for z G C, with 3f(z) > 0, we define the 
matrix-valued functions g&(z) = Pa(H a — z)~ 1 Pa and g&(z) = Pa(Ha — z)~ 1 Pa- 



3.1 Schur 's Complement and Krem's Formula 

Lemma 2 The following formula holds 

1 

QA\Z) = — 7 Krein's formula 

VA + gAiz)- 1 

Proof: The Schur complement formula [18J (also called Feshbach's projection method 
[60) states that if H = H* is a bounded, selfadjoint operator on some Hilbert space and 
if P = 1 — Q is an orthonormal projection, then 



H eff {z) = PHP + PHQ—— QHP . (17) 



H-z H ef f(z)-zP' " w QHQ 



1 The Schur complement method [18j is widely used in numerical analysis under this name, while Math- 
ematical Physicists prefer the reference to Feshbach[6]. It is also called Feshbach-Fano [5] or Feshbach- 
Lowdin [12j in Quantum Chemistry. This method is used in various algorithms in Quantum Chemistry 
(ab initio calculations), in Solid State Physics (the muffin tin approximation, LMTO) as well as in Nuclear 
Physics. The formula used above is found in the original paper of Schur [18j (the formula is on p. 217). 
The formula has been proposed also by an astronomer Tadeusz Banachiewicz in 1937, even though closely 
related results were obtained in 1923 by Hans Boltz and in 1933 by Ralf Rohan [17]. Applied to the Green 
function of a selfadjoint operator with finite rank perturbation, it becomes the Krem formula [llj . 



s 



For completeness we provide a proof of (ITT)) in section [51 Applied to g A (z) = P A (H A — 
z)~ 1 P A gives 

^a(^)- 1 = H A + P A H A P A \ A —^ P A \ A H A P A . 

-Ha\a — z 

By definition, H A = H A +V A while H A \ A = #a\a> so that, applying the Schur complement 
formula to g A (z) instead, gives the desired result 



2i 

1 ^A^)" 1 *-^)-^ 1 



^a + ^a^)- 1 V 2l J V A + g A (z 

' ~ [^a(z)- 1 } ' 



-1 * 



□ 



g A {z)~ l = V A + g A {z) x . 

Lemma 3 If^sz> 0, then ^sg A (z) > and — ^{gA(z)^ 1 } > 0. 
Proof: The resolvent equation gives 



^SZ ^SZ 

%9a(z) = P A — t^Pa > 0, %~g A (z) = P ATT -P A > 0. 

\Ha-z\ \H a -z\ 2 

Using A^ 1 — A^ 1 * = A' 1 {A* — A}A~ 1 * gives the other inequality. □ 

Lemma 4 The following formula holds: 

a +/o ( w det{-^ A (2:) _1 } 
det{% A (^)} = I j f - - 7 W1|2 
\det{V A + g A {z) L }\ 2 

Proof: By definition of the imaginary part, using Lemma [2] gives 



o / x 9a{z) -g&{z 



V A + g A (z)~ 1 1 J " v '' J V A + ~g A {z)~ 1 * 
Taking the determinant of both sides gives the result. □ 

Corollary 1 If Eg denotes the average over the potentials V x for x G S, the following 
estimate holds: 



E A (det{% A (2)}) < E AXA f || /0 ||S o det{-^ A (^)- 1 } / dV A - 



1 



det{V A + g A (z)-^ 



Proof: By definition and since |A] = n, if / is a nonnegative function of V = (V x ) x£ \ 
then 

Mf) = [ \[dV xP {V x )f{V) < \\p\rjZA\A ( [ i[dV x f(V)). 

Since <?a(z) does not depend on Va, the result follows from Lemma HI □ 

Lemma 5 Let M be a complex n x n matrix such that M = B — iA with B = B* and A 
positive definite. Then, taking the principal branch of the square root, 

1 e W4 f e -MMu)/2 (lg) 



v/delM 7 R n (2vr)™/ 

Proof: Since A > 0, it follows that iM has a positive definite real part, so that the 
integral converges and is analytic in M. The formula follows from standard Gaussian 
integrals. □ 

Lemma 6 Let F be an integrable function on R 2 xl 2 . Then the following formula holds 

[ d 2 ud 2 v F(u,v) 6 ( u2 ~^ \ = f d 2 u [ d9F(u,R e u), (19) 
where Rg denotes the rotation of angle 9 in 



5 2 



Proof: Let u be expressed in polar coordinates (r, <ft). The change of variable s = u 2 /2 = 
r 2 /2 gives u = (V^s, <fi) and d 2 u = dsd<fi. In much the same way, v can be expressed as 
(y/2~i,xjj). Thus the integral becomes 



/ 



d 2 u d 2 v F(u, v) S 

xR 2 

oo poo p2tt 

ds / dt 

o Jo JO JO 

oo /-27T /•2tt 

ds 

JO JO 



u 2 — v 2 



2k 



dip F(V2s,(f);V2t 1 ip)5(s-t) 
dip F(V2s,0; V2~s,ip) . 



Setting ip = 9 + (ft gives the result. 



□ 
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3.2 Proof of Theorem H 

Thanks to Corollary [H the main Theorem [2] follows from the following inequality. 
Lemma 7 The following estimate holds 

f 1 7T™ 



/ r a Idetj^A + ^A^)- 1 }! 2 " det{-^ A (^)- 1 } ' 
Proof: Using the Gaussian integral in Lemma [5J the integral J can be written as 

J= [ dV A f ^1^2 d^d^ c _ 1/2 y^ i2 {(uMVA+~9A(z)- 1 )n i )-(vMVA+9A(z)- 1 )v i )} 

V 7 (r a )X 4 (2tt)« (2tt)« 

(20) 

Let = -U2(^)) £ R 2 where Wj = (wi^))^^ e I n mucn the same way let 

v(x) = (vi(x),V2(x)) G M. 2 be used in this integral. The term V x appears in the Gaussian 
exponent with the factor (— t/2)V x (u(x) 2 — v(x) 2 ). Hence integration over V x gives 



/ 



dV x e -^K-) 2 -^) 2 )/2 = 27T5 ( ' 7(X)2 ~' ;(X) ' 



Inserting this result in eq. ( 1201) . using Lemma [6] leads to 

J=TT /" 27r ^. f d<2nU c - 1 / 2 i(»l».9A(^)- 1 g)+(fi(e)g|^A(^- 1 fi(e)») i 
4i^o " y(MA)x2 (2tt)" 

where -R(#) is the orthogonal 2n x 2n matrix acting on u = {u{x)) x£ a by 

(R(6)u)(x) = R 8x u(x). 

Because of Lemma [3l we know that — S^a(-z) 1 > 0, so that the Gaussian term can be 
bounded from above by 

j<TT [ 2 * d6 x [ d2n ^ c -i/2{(g|(-3fl A ( 2 )- 1 )g)+we)g|(-3gAU)- 1 )fl(o)g) _ 

Thus a Schwarz inequality, the rotational invariance of the measure d 2n u and another use 
of the Gaussian formula given in Lemma [5] gives the bound 

7<TT df) f ^ip-ISK-afeW- 1 )*? _ 7r " 
" lih ' V )x2 (2vr)« e ^ det{-^ A (,)-i} ' 

proving the theorem. □ 
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4 Applications of the Correlation Estimate and Re- 
lated Results 



4.1 Level Statistics and Simplicity of Eigenvalues 

The proof of Poisson statistics for the eigenvalue level spacing in the thermodynamic limit 
for general Anderson Hamiltonians Hq + V as in (J9j) follows as in Minami's article provided 
several other conditions are satisfied. In addition to the selfadjointness and boundedness 
of H , we require that H be translation invariant, so that the DOS exists, and that the 
off-diagonal matrix elements of H decay exponentially in \x — y\ with a uniform rate. In 
addition to the positivity of the DOS at energy E, discussed in the introduction, Minami 
requires the exponential decay of the expectation of a fractional moment of the Green's 
function of H u . Let us describe this fractional moment condition. The Green's function 
G\(x,y; z) for the restriction of to a finite cube A C Z d with Dirichlet boundary 
conditions is required to satisfy the following bound. There is some < s < 1 and 
constants C s > and > so that 

E{|G A (a;, y; E + *e)| s } < C s e~ a ^ x - y \ (21) 

provided x G A, y G <9A, and z G {w G C | > 0, \w — E\ < r}, for some r > 0. 

Corollary 2 Consider the general Anderson model |PP with a bounded, translation-invariant, 
selfadjoint H having matrix elements satisfying \ (x\H \y) \ < Ce~ v \ x ~ y \ for some C < oo 
and tj > 0. Suppose that the DOS n(E) for exists at energy E and is positive. Suppose 
also that the expectation of some fractional moment of the Green's function of decays 
exponentially fast as described in (fUJ). Then, the point process (Tj|) converges weakly to 
the Poisson point process with intensity measure n(E) dx. 

We don't provide a detailed proof of this here, as it is easily checked that under the as- 
sumption of exponential decay of |(a;|iio|2/)| the remaining arguments of Minami's proof 
of Poisson statistics go through. Translation invariance of H comes in as an extra as- 
sumption to guarantee ergodicity of H u , and thus existence of the IDS (j2J). 

Exponential decay of | (x\H \y) \ implies the strong localization condition (j2"T]) at ex- 
treme energies or high disorder [2], or, for low disorder, at band edges pQ. Also, exponen- 
tial bounds of the form ( 1211) imply almost sure pure point spectrum for the energies at 
which they hold and exponential decay of the corresponding eigenfunctions. 

These conditions on Hq and the decay estimate (1211) also insure that the result of 
Klein and Molchanov [10] (which uses [2] and thus rapid off-diagonal decay of the matrix 
elements of H ) on the almost sure simplicity of eigenvalues are applicable in the above 
situation, thus: 

Corollary 3 The eigenvalues of the general Anderson model considered in Corollary [H 
in the region of localization are simple almost surely. 
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In fact, all of the above can be extended to Hq with sufficiently rapid power decay of 
the off-diagonal elements. The works [2] and pQ discuss how a result somewhat weaker 
than (12T|) can be obtained in this case. In particular, this only gives power decay of 
eigenfunctions, which for sufficiently fast power decay still allows to apply the result of 
[TO] . Moreover, a thorough analysis of Minami's proof shows that it works for suitable 
power decay. 



4.2 A different proof of Theorem 2 

After we finished the proof of Theorem [2J we received the preprint of Graf and Vaghi 
[7] in which they proved essentially the same result using a different approach. One of 
their main motivations was to eliminate Minami's symmetry condition on the Green's 
function Ga(x,U]z) = Ga(v,x;z) thus allowing magnetic fields as in ([Til . They base 
their calculation on the following lemma. By diag (v±, . . . , v n ), we mean the n x n-matrix 
with only nonzero diagonal entries V\, . . . , v n . 

Lemma 8 Let A = (a^) be an n x n matrix with > 0. then, 

J dv x --- dv n det(9f[ diag (v u ...,v n )- A]' 1 ) < rr n . (22) 

It is not surprising that the proof of Lemma [8] involves the Schur complement formula. 
One applies Lemma [8] by noting that the argument of the determinant on the left side of 
(prj (for the case n = 2) may be written as the imaginary part of a matrix of the form 
[ diag (v\, . . . ,v n ) — A]~ l by Krein's formula, where A is obtained from H\ by setting 
V(x) = V(y) = 0. Without explicitly stating it, Graf and Vaghi also indicate that a 
bound as in f|T2l follows from (1221) for general n. 

The key to proving Lemma [8] for n = 2 are the integral formulas 

1 71 

dx -, -777 = = , (23) 

\ax + b\ 2 Q(ba) 
assuming a, b G C and Q(ba) > 0, and 

dx — — = - (24) 

ax z + bx + c V4ac - b 2 

assuming a, b, c e K, a > 0, and 4ac — b 2 > 0. The case for general n is obtained by 
induction. 



4.3 Joint Energy-Space Distributions 

We mention two related results of interest. Nakano [15] recently obtained some quanti- 
tative results providing insight into the idea, going back to Mott, that when eigenvalues 
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in the localization regime are close together, the centers of localization are far apart. 
Roughly, Nakano proves that for any subinterval J of energies in the localization regime 
with sufficiently small length, there is at most one eigenvalue of H w in J with a localiza- 
tion center in a sufficiently large cube about any point with probability one. His proof 
uses Minami's estimate in the form ([6l) and the multiscale analysis. In this sense, the 
centers of localization are repulsive. On the other hand, if one studies an appropriately 
scaled space and energy distribution of the eigenfunctions in the localization regime in the 
thermodynamic limit, Killip and Nakano [9] proved that this distribution is Poissonian, 
extending Minami's work for the Anderson model ([1]). They define a measure d£ on R d+1 
by the following functional. For / G Cq(R) and g G Co(R d ), consider the map 

/, g - tr(f(H)g(-)) = [ f(E)g(x) d£(E, x). (25) 



This measure is supported on E x Z d , where E C 1 is the deterministic spectrum of H u . 
They perform a microscopic rescaling of d^ in both energy and space to obtain a measure 
d^i as follows 

J f(E,x) d£ L (E,x) = J f{L\E-Eo),xL- x ) d£(E,x), (26) 

where E is a fixed energy for which the density of states n exists and is positive. In the 
limit L —>■ oo, they prove that this rescaled measure converges in distribution to a Poisson 
point process on R x R d with intensity given by n(Eo)dE x dx. This work also relies on 
Minami's estimate (jSJ) but uses the fractional moment estimates rather than multiscale 
analysis. Both of these papers treat the standard Anderson model (pQ) so that Theorem 2 
extends the results to more general lattice operators of the form H + V^. 



5 Appendix: The Schur Complement Formula 



We prove the Schur complement formula for a selfadjoint operator H and an orthogonal 
projection P with Q = 1 — P on a Hilbert space Ti. In the case that H is unbounded, we 
assume that PH. C D(H). Let zeC and suppose that Q(H — z)Q is boundedly invertible 
on the range of Q (as always the case for z G C \ R). We write Rq{z) = (Q(H — z)Q)~ x 
for the resolvent of the reduced operator. We write H — z on Ti as the matrix 



H-z = 

We introduce the triangular matrix L given by 



P(H-z)P PHQ 
QHP Q(H - z)Q 



(27) 



P 

-R Q {z)QHP R q {z) 



(2* 
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The Schur complement of Q(H — z)Q is defined as S{z) = P(H — z)P — PHQRq(z)QHP. 
We assume that S(z) is boundedly invertible on the range of P (true for z G C \ M.). 
Multiplying (H — z) on the right by L we obtain 



(H-z)-L 



S(z) PHQR Q (z) 
QP(H-z)P 



(29) 



This matrix may be inverted, and multiplying the inverse ((H — z)L) 1 = L 1 R(z) on 
the left by L gives 



R(z) 



S(z)- 1 -S{z)- l PHQR Q {z) 
-R Q {z)QHPS{zY 1 Rq(z) + R Q (z)QHPS(z)- 1 QHPR Q (z) 



(30) 



The formula for PR(z)P readily follows from (I3"0"|) since in matrix notation P is block 
diagonal. 
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